Kinetic exchange models: From molecular physics to social science 
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We discuss several multi-agent models that have their origin in the kinetic exchange theory of 
statistical mechanics and have been recently applied to a variety of problems in the social sciences. 
This class of models can be easily adapted for simulations in areas other than physics, such as the 
modeling of income and wealth distributions in economics and opinion dynamics in sociology. 



I. INTRODUCTION 

The application of probability theory to the study of 
gases led to the formulation of the kinetic theory of gases, 
prepared the basis for the formulation of the Maxwell ve- 
locity distribution, and later, for the development of sta- 
tistical mechanics. ■It'S. The initial triumph was that the 
empirical laws of thermodynamics were obtained from 
the basic assumptions of statistical mechanics. Today 
we understand that the scope of statistical mechanics is 
much broader. Because it has the tools to treat macro- 
scopic systems with a large number of microscopic con- 
stituents, it can be naturally applied to many types of 
systems. 

Kinetic exchange models, the subject of this article, are 
stochastic models which have a straightforward interpre- 
tation in terms of energy exchanges in a gas. However, 
they can be suitably adapted and used to study problems 
in the social sciences, as illustrated by recent work.— 
The applications of kinetic exchange models to problems 
in fields from molecular physics to the social sciences il- 
lustrates the historical link between statistical mechanics 
and the social sciences. Recently, this link has been redis- 
covered due to the development of complex systems the- 
ory and social dynamics. In fact, statistics, which was an 
important basis for Maxwell's and Boltzmann's work and 
the foundation of statistical mechanics, originated from 
the study of demographic datai^ The idea that a large 
number of molecules and a large social group have many 
important common features, in particular, that they are 
predictable due to the high number of their components 
despite their intrinsic random character, was shared by 
many investigators. Boltzmann wrote that "molecules 
are like so many individuals, having the most various 
states of motion," when writing about the foundations of 
statistical mechanicsi^ 

Multi-agent models are a class of models where the ac- 
tions and interactions of autonomous agents, which may 
represent individuals, organizations, societies, etc., can 
be used to understand the behavior of the system as a 
whole. The simple formulation and numerical implemen- 
tation of kinetic exchange models make them attractive 
in many disciplines. They can be regarded as minimal 



prototypical models of complex systems consisting of a 
set of (possibly) heterogeneous units interacting accord- 
ing to simple laws, yet they are able to exhibit well de- 
fined states with robust probability distributions. 

In this article we focus on the formulation, interpre- 
tation, and simulation of a few representative examples 
of kinetic exchange models. Other social science applica- 
tions of kinetic exchange models are also summarized. 



II. HOMOGENEOUS KINETIC EXCHANGE 
MODELS 



Kinetic exchange models in molecular physics: 
Thermal relaxation in d-dimensions 



We introduce the general structure of a kinetic ex- 
change model by a simple example. It is assumed that the 
N (minimally) interacting units {i}, with i — 1,2, . . . , N , 
are molecules of a gas with no interaction energy and the 
variables {wi} represent their kinetic energies, such that 
Wi > 0. The time evolution of the system proceeds by 
a discrete stochastic dynamics. A series of updates of 
the kinetic energies Wi(t) are made at the discrete times 
t = 0,1, . . . . Each update takes into account the effect of 
a collision between two molecules. The time step, which 
can be set to At = 1 without loss of generality, repre- 
sents the average time interval between two consecutive 
molecular collisions; that is, on average, after each time 
step At, two molecules i and j undergo a scattering pro- 
cess and an update of their kinetic energies Wi and Wj is 
made. 

The evolution of the system is accomplished by the 
following steps at each time t. 

1. Randomly choose a pair of molecules i and j {i ^ j 
and 1 < i,j < N), with kinetic energies Wi and wj, 
respectively; they represent the molecules undergo- 
ing a collision. 

2. Compute the net amount Awij of kinetic energy to 
be exchanged between the molecules; its value is a 
function of the initial kinetic energies Wi and Wj 
and depends on the model considered. 
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3. Perform the energy exchange between i and j by 
updating their kinetic energies, 
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The total kinetic energy is conserved during an in- 
teraction. 

4. Set t ^ t + 1 and go to step 1. 

The form of the function Awij defines the specific 
model. Kinetic exchange models describe the dynamics 
at a microscopic level, based on single molecular colli- 
sions. Such a representation can be optimal in terms of 
simplicity and computational efficiency when the focus is 
on the energy dynamics, because particles are described 
by their energy degree of freedom w only, rather than by 
the entire set of their 2d position and momentum coor- 
dinates, for a d-dimensional system. 

As our first simple example, consider the reshuffling 
rule: Wi — >■ e{wi + Wj), Wj — >■ (1 — e){wi + Wj), where 
e is a stochastic variable drawn as a uniform random 
number between and 1. This rule corresponds to 
Awij = (1 — e)wi — ewj. In this case, the algorithm we 
have outlined leads from arbitrary initial conditions to 
the Boltzmann-Gibbs energy distribution at equilibrium, 
/(w) = /3exp(— /3w), where /? = l/(w) and (w) repre- 
sents the mean energy of a single molecule. The the- 
oretical derivations of this result using the Boltzmann 
transport equation, or entropy maximization principle, 
or simple probabilistic arguments, can be found in most 
standard textbooks of statistical mechanics. 

As a more general example, consider the relaxation in 
energy space of a gas in d-dimensions. We assume that 
d > 1 because the momentum and energy distributions 
of a one-dimensional gas (where only head-on collisions 
occur) do not change with time. Although the model can 
be conceived most easily for a gas in three dimensions, 
the most surprising features of kinetic exchange models 
appear for a gas in an arbitrary number of dimensions, a 
case relevant for the treatment of heterogeneous systems 
in Sec. IIIII For a gas in d-dimensions the form of the 
update rule, and in particular of Awij, can be derived 
exactly from energy and momentum conservation during 
a collision between two particles i and j. If the respective 
d-dimensional vectors of the particle initial momenta are 
Pi and Pj, we find^'' 
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where cos ak is the direction cosine of momentum 
{k = with respect to the direction of the transferred 
momentum Ap^ = Pi — Pj- The directions of the two 
colliding particles can be assumed to be random using 
the hypothesis of molecular chaos 



We can now study the time evolution by randomly 
choosing at each time step two new values for rj, in 
Eq. ([2]), instead of maintaining a list of momentum co- 
ordinates, as is commonly done in a molecular dynamics 
simulation. Then we use Eq. ([l]) to compute the new 
particle energies Wi and wj. Note that the r^'s are not 
uniformly distributed in (0, 1), and thus the form of their 
probability distribution function has to be chosen with 
some prudence. In fact, their distribution strongly de- 
pends on the spatial dimension d, their average value 
being (r^) = 1/d; see Ref. [13 for further details. This 
dependence of (rk) on d can be intuitively understood 
from kinetic theory: the greater the value of d, the 
more unlikely it becomes that rk assumes values close 
to Tk = 1 (corresponding to a one-dimensional-like head 
on collision). Hence, a simple and computationally effi- 
cient choice is a uniform random distribution f{rk), lim- 
ited in the interval (0,2/d), such that the average value 
{rk) = 1/d. 

Simulations of this model system using random num- 
bers in place of the ri's in Eq. for d = 2, give 
the equilibrium Boltzmann-Gibbs distribution: f{w) = 
/3exp(— /3u'), where /? = l/(w), as before. For d > 2, 
we obtain the d-dimensional generalization of the stan- 
dard Boltzmann distributionjiSr— namely the Gamma 
(r) distribution^"^'^^ characterized by a shape parameter 
a equal to half of the spatial dimension. 



f{w,a,9) = 



d/2 

(w) I a . 



(5) 

(6) 
(7) 



The scale parameter d of the F-distribution is fixed, by 
definition, by Eq. ([7])J^ii^ From the equipartition the- 
orem in classical statistical mechanics, w = dk^T/2. 
Hence, we see that Eq. ([7]) identifies the scale parameter 
9 as the absolute temperature (in energy units) given by 
9 = k-oT —1/(3. Therefore, the same Boltzmann factor, 
exp(— i/7/0), is present in the equilibrium distribution in- 
dependently of the dimension d, and the prefactor w°'~^ 
depends on d, because it takes into account the phase- 
space volume proportional to p'^ oc w'^^'^, where p is the 
momentum modulus. 



B. Modeling the wealth distribution 

Besides their interpretation in the context of physics, 
kinetic exchange models are applicable to problems in the 
social sciences such as opinion dynamics and the model- 
ing of economic systems. The latter models demonstrate 
their adaptability and promise, particularly in their ap- 
plication to wealth distributions. In this section we con- 
sider how some kinetic exchange models describe the for- 
mation of the wealth distribution observed in a society, 
as a consequence of binary wealth exchanges between in- 
dividuals. We assume here that the standard definition 
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FIG. 1. Cumulative probability distribution of the net 
wealth, composed of assets (including cash, stocks, property, 
and household goods) and liabilities (including mortgages 
and other debts) in the United Kingdom shown on log-log 
(main panel) and log-linear (inset) scales. Points represent 
the data from the Inland Revenue, and solid lines are fits 
to the Boltzmann-Gibbs (exponential) and Pareto (power) 
distributions.— 

of wealth is the set of all those things with some mone- 
tary or exchange value. The shape of a typical wealth 
distribution is complex, with a Boltzmann-Gibbs (ex- 
ponential) behavior at lower and intermediate values of 
wealth, w < Wc, and a power law (Pareto) tail at the 
larger wealth values, w > Wc, where Wc is a crossover 
value that depends on the numerical fitting of the data. 
The Pareto law is expressed as: 

f{w) ^ {w > Wc), (8) 

where f{w) is the probability density, and the exponent a 
is the Pareto exponent, which has a value between 1 and 
3. Reference |l6| gives concise accounts of empirical, nu- 
merical and analytical studies on this subject. Figure [T] 
shows a plot of the cumulative probability distribution 
for wealth data. It is more practical to plot the cumu- 
lative probability distribution C{w) = f{w') dw' as a 
function of the wealth w, rather than f(w), because sta- 
tistical data is usually reported at non-uniform intervals 
of w. Interestingly, when /(w) is an exponential or a 
power-law function, then the respective C{'w) is also an 
exponential or a power-law function. For more details of 
kinetic wealth exchange models, see Refs.Si, andS 

Among the first examples of kinetic exchange models 
of markets proposed by non-physicists, we mention the 
work of Angle, a sociologist. ^^"^'^ In the first version of the 
modeliiii^ the wealth exchanged at a encounter between 
two agents has the form of Eq. ([2]) , where = tnr] and 
rj = (1 — e) nrj. Here k represents the maximum fraction 
of wealth that can be lost, e is a uniform random number 
in (0, 1), and ?] is a stochastic variable with the values 



or 1 with a probability distribution depending on the 
difference between the wealth of the two agents Wi — Wj . 
In the second model by Angle, referred to as the one- 
parameter inequality procesSfi^iSS n = kt] and rj = k(1 — 
rj), with a similar meaning for k and rj. In contrast, in 
the basic version of the models^ir— introduced earlier by 
the economist, Bennati, Awij is a constant, independent 
of other parameters. 

In the models introduced by physicists, we first con- 
sider the work of Dragulescu and YakovenkOj-^^ where 
ri — 1 — rj e, corresponding to a random reshuffling 
scheme of the total wealth of the two agents. We next 
consider a simple prototypical kinetic exchange model 
with a F-function equilibrium distribution, which was 
introduced by Chakraborti and Chakrabarti^^ to de- 
scribe the trade activity between N entities (for example, 
agents, firms, and companies) employing a saving crite- 
rion when carrying out their trades. Trades are repre- 
sented by pair-wise wealth exchanges with the dynamics 
given by Eq. ([ij with 

Awij = K{lwi - ewj) = (1 - A)(ewi - ewj) , (9) 

where e = 1 — e. The exchange parameter k £ (0, 1) (or 
the saving parameter A = 1 — k) defines the maximum 
fraction of the wealth w in the exchange process (or the 
minimum fraction of w preserved during the exchange). 
The parameter k (or A) also determines the time scale of 
the relaxation process as well as the mean value (w) in 
equilibrium;^ The equilibrium wealth distribution of the 
system is well described by the F-distribution in Eq. ([S|) , 
where a fit suggests the empirical formula 

1 + 2A 3 
1 — A K 

which relates the shape parameter a to the saving pa- 
rameter A. Note that for A (or k 1), a — >■ 1, 
corresponding to the exponential function. In all the 
various models we have mentioned, the equilibrium dis- 
tribution is well fitted by a F-distribution. This shape of 
the equilibrium distribution shows good agreement with 
empirical wealth distributions at small and intermediate 
valuesi^i^ (see Fig. [1]). 

Numerical results of the model defined by Eq. (jH]) are 
compared in Fig. [2] with a fit based on the F-distribution. 
They were obtained using 10'' agents and lO'' time steps. 
The convergence to equilibrium is fast and such a long 
simulation time was used only to accumulate statistics 
(every 100 time steps). Note that here one time step 
is actually a loop over N exchanges. The number of 
agents employed has to be chosen sufficiently large to 
ensure enough statistics — implying a good quality of 
the wealth histogram — at the desired frequency scale, 
otherwise a very irregular histogram may appear at the 
smallest and largest values of w. The value N = 10^ is 
still manageable computationally, yet sufficiently large to 
clearly see frequencies as small as 10~^. 
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FIG. 2. (Color online) Equilibrium wealth distributions of 
the model introduced in Ref. [2^ and defined by Eq. (|9} for 
different values of the saving parameter A for linear (top) and 
semi-log (bottom) scales. 



C. Other social science applications 

Lallouache et aL - constructed a homogeneous model 
for the collective dynamics of opinion formation in a so- 
ciety by modifying the kinetic exchange dynamics stud- 
ied in the context of markets. In this model the opinion 
of the ith agent is represented by a continuous variable 
Wi € (—1,-1-1). Interactions between agents who modify 
their opinions are assumed to take place through random 
two-body encounters, with a dynamics defined by 



X[wi + ewj] 
X[wj + e'wi] 



(11a) 
(lib) 



where A € (0, 1) is the (constant) "conviction" parameter 
and e, e' £ (0, 1) are uncorrelated uniformly distributed 
stochastic processes. Note that there is no step-wise 
conservation of opinion, unlike for the preceding wealth 
models. Remarkably, it is found that there is an ap- 
pearance of polarity or consensus, starting from initial 
random disorder (where the Wi are uniformly distributed 
with positive and negative values). In the language of 
physics, there is a "spontaneous symmetry-breaking tran- 
sition" in the system: starting from the mean value of Wi , 



iZ'(O) = 0, the system evolves either to the "para" state 
with w{t > r) = 0, where all agents have the opinion 
zero (for values of A < 2/3), or to the "symmetry bro- 
ken" state, with w{t > r) ^ 0, where all the agents have 
either all positive or all negative opinions (for A > 2/3). 
The time t > r, where t is the relaxation time for the 
system. The relaxation behavior of the system shows a 
critical divergence of r at A = Ac = 2/3. Sen^* gener- 
alized the model by Lallouache et al. by introducing an 
additional parameter /x to represent the influencing abil- 
ity of individuals and studied the corresponding phase 
transitions. 

In a different contextj_similar in spirit to the work pre- 
sented in Refs. and 113, Ghosh et al.'^- considered an 
economic model in which there is a poverty threshold 
6 > such that at any time t, at least one of the two 
interacting agents is "poor;" that is, its wealth satisfies 
w < 9. The central role assigned to the poor traders 
produces various new features such as a different form of 
the equilibrium distribution and a phase transition in the 
fraction of poor agents as a function of 9. These models 
of kinetic exchanges, which draw inspiration from vari- 
ous socio-economic contexts, may also serve the purpose 
of introducing the ideas of phase transitions and critical 
phenomena in statistical physics. 



III. HETEROGENEOUS KINETIC EXCHANGE 
MODELS 

An interesting generalization of the homogeneous ki- 
netic exchange models we have discussed so far is the in- 
troduction of heterogeneity. Probably the most relevant 
applications of heterogeneous kinetic exchange models in 
the social sciences is the prediction of a realistic shape for 
the wealth distribution which includes the Pareto power 
law at the largest wealth values (see Fig.[T]and Sec. Ill Bl) . 

We consider again the model defined by Eq. ^ and in- 
troduce heterogeneity by diversifying the parameter k or 
equivalently the saving parameter A, meaning that each 
term KWi (or Xwi) is replaced by KiWi (or XiWi), thus 
obtaining'^— 

Awij ~ e KiWi ^ eKjWj = — Xi)wi — e{l^ Xj)wj . (12) 

As a simple example, we consider a set of heteroge- 
neous agents with parameters Ki uniformly distributed 
in the interval (0,1). By repeating the simulations us- 
ing Eq. ([T2|) . it is found that the shape of the separate 
equilibrium wealth distributions fi{'w) of each agent is 
still a F-distribution. However, there is a surprise in 
the wealth distribution of the system f{w), given by 
the sum of the wealth distributions of the single agents, 
fiw) = "^ifiiw)- As other analytical and numerical 
studies have also shown, f{w) has an exponential form 
until intermediate w- values, and a Pareto power law de- 
velops at the largest values of w (see Fig. [3l).'^-i^ Such a 
shape is similar to real wealth distributions such as that 
shown in Fig. [TJ This shape of the equilibrium wealth 
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FIG. 3. (Color online) Wealth distribution f{w) for uniformly 
distributed (or Ai) in the interval (0,1); f{w) is decomposed 
into partial distributions /i(w), where each /i(w) is obtained 
by counting the statistics of those agents with parameter Xi 
in a specific sub-interval (from Ref. iSfj). (a) Decomposition 
of f{w) into ten partial distributions in the A-subintervals 
(0, 0.1), (0.1, 0.2) ...(0.9, 1). (b) The last distribution of 
(a) in the A-interval (0.9, 1) is decomposed into partial dis- 
tributions obtained by counting the statistics of agents with 
A-subintervals (0.9, 0.91), (0.91, 0.92) . . . (0.99, 1). Note how 
the power law appears as a consequence of the superposition 
of the partial distributions. 



distribution f{w) is robust with respect to the details of 
the system and the values of the other parameters, as 
long as the values of the Hi are sufficiently spread over 
the entire interval k = (0, 1). It is the group of agents 
with K w (A « 1) that are crucial for the appearance 
of a power law. Not all agents have to differ from each 
other, as is best illustrated by repeating the simulation 
using a different distribution for the K-parameters or the 
As, in which 99% of the agent population has A — 0.2, 
and only 1% of the population is heterogeneous with A 
in the interval (0, 1) (see Fig. S]). 

The heterogeneous model necessarily uses a finite up- 



w 

FIG. 4. (Color online) Example of a realistic wealth distri- 
bution, from Ref. |35| . Continuous curve: wealth distribu- 
tion obtained by simulations of a mixed population of agents, 
such that 1% of the agents have uniformly distributed sav- 
ing propensities Xi £ (0, 1) and the other 99% of the agents 
have Ai = 0.2. Dotted curve: exponential wealth distribution 
oc exp{—w/{w)), with the same average wealth, plotted for 
comparison with the distribution in the intermediate-wealth 
region. Dashed curve: Pareto power law cx plotted for 

comparison with the large-income part of the distribution. 



per cutoff Amax < 1, when considering the saving pa- 
rameter distribution, which directly determines the cutoff 
Wmax of the wealth distribution, analogous to the cutoff 
observed in real distributions: the closer Amax is to one, 
the larger Wmax and the wider the interval in which the 
power law is observed. "^^ 

The A-cutoff is closely related to the relaxation process, 
whose time scales for agent i is proportional to 1/(1 — 
Aj)3£ Thus the slowest convergence rate is determined 
by 1 — Amax- The finite A-cutoff used in simulations of 
heterogeneous kinetic exchange models is not a limitation 
of the model, but reflects an important feature of real 
wealth distributions. 



IV. COMMENTS 

There are several analytical studies which complement 
the simulations and empirical studies. We refer the 
reader to Ref. [l^, which gives concise accounts of studies 
on this subject, along with an exhaustive list of refer- 
ences. 

The dynamics of kinetic exchange models are some- 
times criticized for being based on an approach that is 
far from an actual economics perspective. However, such 
a dynamics can also be derived from microeconomics 
theoryj22, Although standard economics theory assumes 
that the activities of individual agents are driven by 
the utility maximization principle, the alternative picture 
that we have described is that the agents can be viewed 
as particles exchanging "wealth," instead of energy, and 
trading in wealth (energy) conserving two-body scatter- 
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ing, as in entropy maximization based on the kinetic the- 
ory of gases. -'^ This qualitative analogy between the two 
maximization principles is not new - both economists 
and physicists had noted it in many contexts, but this 
equivalence has gained firmer ground only recentlyi^ 

V. SUGGESTED PROBLEMS 

The wealth exchange models, we have discussed are 
minimal and lend themselves to various generalizations 
and modifications to study the effect of additional fea- 
tures. Choose the kinetic exchange model, for instance, 
defined by: Wi ^ Wi — Awij and wj wj + Awij , with 
Awij = (1 — e)wi — ewj. Perform the simulations and 
check whether the corresponding equilibrium Boltzmann- 
Gibbs distribution f{w) — (3 exp{— f3w) , where (3 — 1/ (w) 
is reached, starting from any arbitrary initial distribution 
of w amongst the agents. Next, perform the simulations 
for Wi —7- Wi — Awij and Wj — > Wj + Awij, when the 
form of Awij or the dynamical equations are modified, 
as suggested in the following. 

• Different exchange rules: Verify that the 
final equilibrium wealth distribution remains 
Boltzmann-Gibbs distribution, when employing the 
rule that at each step an arbitrary constant amount 
[Awij = wq] is exchanged, or at each step a fraction 
of the average wealth in the system is exchanged)^ 
However, if instead the minimum wealth between 
the two agents is used as the exchanged amount 
[Awjj = mhi{wi,Wj)], condensation of the total 
wealth in the hands of a single agent will take 
place 1^ 

• Taxation on transactions: Modiiy the dynami- 
cal equations by adding a term that takes as tax. 



a fixed amount of wealth Swq or a fixed percent- 
age of the wealth exchanged Awij from agents i 
and j , and redistributes it uniformly among all the 
others i2i 

• Debt: Modify the dynamical equations such that 
the agents are allowed to borrow or take loans up to 
a maximum amount Wmax so that their wealth may 
become negative but not smaller than — u;,nax- This 
modification can be implemented by increasing the 
total wealth by an additional amount Wmax"^ 

• Fixed saving: Add to the model a fixed saving wq 
for each of the two agents, so that the total wealth 
is diminished by an amount 2wo^ 

• Many-agent interactions: Set up a new dynami- 
cal rule in which wealth is redistributed in encoun- 
ters between more than two agents, for example, 
among three agents z, j, and k. The other de- 
tails of the transaction should be similar to that 
of a pair- wise transaction between agents i and j. 
Some care has to be taken to properly define the 
respective random fractions of wealth assigned to 
each agent to make the process symmetrical with 
respect to all the agents involved in the encounter. 
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